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Abstract 

This article investigates the predictions of an inflationary phase starting from a homogeneous 
and anisotropic universe of the Bianchi / type. After discussing the evolution of the background 
spacetime, focusing on the number of e-folds and the isotropization, we solve the perturbation 
equations and predict the power spectra of the curvature perturbations and gravity waves at the 
end of inflation. 

The main features of the early anisotropic phase is (1) a dependence of the spectra on the 
direction of the modes, (2) a coupling between curvature perturbations and gravity waves, and (3) 
the fact that the two gravity waves polarisations do not share the same spectrum on large scales. 
All these effects are signihcant only on large scales and die out on small scales where isotropy is 
recovered. They depend on a characteristic scale that can, but a priori must not, be tuned to some 
observable scale. 

To hx the initial conditions, we propose a procedure that generalises the one standardly used in 
inflation but that takes into account the fact that the WKB regime is violated at early times when 
the shear dominates. We stress that there exist modes that do not satisfy the WKB condition 
during the shear-dominated regime and for which the amplitude at the end of inflation depends on 
unknown initial conditions. On such scales, inflation loses its predictability. 

This study paves the way to the determination of the cosmological signature of a primordial 
shear, whatever the Bianchi / spacetime. It thus stresses the importance of the WKB regime to 
draw inflationary predictions and demonstrates that when the number of e-folds is large enough, 
the predictions converge toward those of inflation in a Friedmann-Lemaitre spacetime but that 
they are less robust in the case of an inflationary era with a small number of e-folds. 
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I. INTRODUCTION 

Inflation [H, [2| (see Ref. for a recent review of its status) is now one of the corner- 
stones of the standard cosmological model. In its simplest form, inflation has very definite 
predictions: the existence of adiabatic initial scalar perturbations and gravitational waves, 
both with Gaussian statistics and an almost scale invariant power spectrum 0, HI- Various 
extensions, which in general involve more fields, allow e.g. for isocurvature perturbations j^, 
non-Gaussianity 0], and modulated fluctuations All these features let us hope that fu- 
ture data will shed some light on the details (and physics) of this primordial phase of the 
universe. 

Almost the entire literature on inflation assumes that the universe is homogeneous 
and isotropic while homogeneity and isotropy are what inflation is supposed to explain. 
Indeed, the dynamics of anisotropic inflationary universes has been widely discussed jof. 
It was demonstrated that under a large variety of conditions, inflation occurs even if the 
spacetime is initially anisotropic [13], regardless of whether it is dominated by a pure 
cosmological constant or a slow-rolling scalar field. The isotropization of the universe was 



even generalized to Bianchi braneworld models recently. It should be emphasized, 
however, that a deviation from isotropy [lo| or flatness [l2| may have a strong effect on the 
dynamics of inflation, and in particular on the number of e-folds. 

The study of the perturbations during the isotropization phase has been overlooked, 
mainly because in the past decades one was mostly focused on large field inflationary models, 
which generically give very long inflationary phases. This was also backed up by the ideas of 
chaotic inflation and eternal inflation [3| . In such cases, it is thus an excellent approximation 
to describe the universe by a Friedmann-Lemaitre (FL) spacetime when focusing on the 
modes observable today since they exited the Hubble radius approximatly in the last 60 e- 
folds. In this case, the origin of the density perturbations is understood as the amplification 
of vacuum quantum fluctuations of the inflaton. In particular, the degrees of freedom that 
should be quantized deep in the inflationary phase, and known as the Mukhanov- Sasaki 



variables [13|, were identified [5|, which completely fixes the initial conditions and makes 
inflation a very predictive theory. 

In the context of string theory, constructing a string compactification whose low energy 
effective Lagrangian is able to produce inflation is challenging (see Ref. [3l)- In particular, 
it has proved to be difficult to build large field models [l5| (in which the infiaton moves 
over large distances compared to the Planck scale in field space). This has led to the idea 
that, in this framework, an infiationary phase with a small number of e-folds is favored 
(see however Ref. [l6|). If so, the predictions of infiation are expected to be sensitive 
to the initial conditions, and in particular the classical inhomogeneities are not expected 
to be exponentially suppressed, which makes the search for large scale deviations from 
homogeneity and isotropy much more motivated, as well as the possibility that the infiaton 
has not reached the infiationary attractor. More important, it is far from obvious (as we 
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shall discuss in details later) that all observable modes can be assumed to be in their Bunch- 
Davies vacuum initially, and more puzzling, that for a given mode modulus the possibility 
of setting the initial conditions will depend on its direction. If so, then the initial conditions 
for these modes would have to be set in the stringy phase, an open issue at the time. 

The theory of cosmological perturbations in a Bianchi universe was roughed out in 
Ref. [13, [H, [3] (see also Ref. [13] and Ref. j2l| for the case of higher-dimensional 
Kaluza-Klein models and Ref. [12] for the quantization of test fields and particle production 
in an anisotropic spacetime). Recently, we performed a full analysis of the cosmological 
perturbations in an arbitrary Bianchi / universe jl^ . It was soon followed by an analysis 
that focused on Bianchi / universe with a planar symmetry. As we shall see in this work, 
the case of a planar symmetric spacetime is not generic (both for the dynamics of the 
background and the evolution of the perturbations). 

From a more observationally oriented perspective, the primordial anisotropy can imprint 
a preferred direction in the primordial power spectra. This could be related to the 
possible large scale statistical anomalies [25| of the cosmic microwave background (CMB) 
anisotropics. Many possible explanations have been proposed, including foregrounds (26j . 
non-trivial spatial topology (23] (which implies a violation of global isotropy (28|), the 
breakdown of local isotropy due to multiple scalar fields [loj , the presence of spinors 3^ or 
dynamical vectors |3l| . the effect of the spatial gradient of the infiaton 3^ or a late time 
violation of isotropy [33|. In the case of universes with planar symmetry, the signatures on 



the CMB were derived by many authors [2J, |3J, |35|. In this case, a large primordial shear 



is indeed necessary, which is not in contradiction with the constraints obtained from the 



CMB [36| or from the big-bang nucleosynthesis [37|. This brings a secondary motivation to 
our analysis: can the CMB anisotropy be related to an anisotropic primordial phase or, on 
the other hand, can it constrain the primordial shear and what are the exact predictions of 
a primordial anisotropic phase? 

This is however not our primary motivation. In the first place, we are interested in 
understanding the genericity of the predictions of inflation, and in particular with respect 
to the symmetries of the background spacetime. As we shall see, the simple extension 
considered in this article drives a lot of questions, concerning both the initial conditions in 
inflation and, more generally, quantum field theory in curved spacetime. 



In this article, we build on our previous work [23| to investigate the dynamics and pre- 
dictions of one field inflation starting from a generic Bianchi I universe. After discussing 
the dynamics of the background in § |Tll we focus on the evolution of the perturbations in 
§ mil In particular, we will need to understand the procedure for the quantization during 
inflation. As we shall see in § llVt this procedure deviates from the one standardly used in 
a Friedmann-Lemaitre spacetime, and even in a planar symmetric spacetime as discussed 



in Ref. [2J]. The main reason for such an extension lies in the fact that there always exist 
modes that were not in a WKB regime during the shear-dominated inflationary era. It 
follows that, while our procedure leads to similar predictions to the standard one on small 
scales, it appears that there is a lack of predictibility on large scales. Indeed, we do not 
want to push our description beyond the Planck or string scale where extensions of general 
relativity have to be considered. This may give a description of both the early phase of the 
inflationary era and also a procedure to fix the initial conditions without any ambiguity. 
In § El we explicitly compute the primordial spectra, for both scalar modes and gravity 
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waves. We will describe in details the effect of the anisotropy and show how the isotropic 
predictions are recovered on small scales. 



II. BACKGROUND DYNAMICS 

We first set our notations in § 111 Al and describe the dynamics of the background, focusing 
on its general solutions in § IIIBI We then turn to the slow-roll regime and to the case of a 
massive free field, that we shall use as our working example, respectively in § 111 CI and § IIIDI 



Definitions and notations 



Bianchi spacetimes [38|] are spatially homogeneous and those of type I have Euclidean 
hypersurfaces of homogeneity. In comoving coordinates, and using cosmic time, their metric 
takes the general form 

3 

ds' = -dt' + J2xKt) {dxf , (2.1) 
1=1 

which includes three a priori different scale factors. It includes the Friedmann-Lemaitre 
spacetimes as a subcase when the three scale factors are equal, and the extensively studied 
planar symmetric universes when only two of the three scale factors are different. The 
average scale factor, defined by 

S{t) ^ [X^{t)X,{t)X,{t)f' , (2.2) 

characterises the volume expansion of the universe. The metric (12. ip can then be recast 
under the equivalent form 

ds^ = -dt^ + S\t)-fij{t)dx'dx^ , (2.3) 

where the "spatial metric", 7jj, is the metric on constant time hypersurfaces. It can be 
decomposed as 

7,,- = exp [2A(t)] Sij , (2.4) 
where the functions /3j must satisfy the constraint 



I]A = 0. (2.5) 



i=l 

Let us introduce some useful definitions. First, we consider the scale factors 

ai = e^'(*) , Xi = Sai . (2.6) 
They are associated to the following Hubble parameters 

H^^, h = ^, A = ^, (2.7) 

O A,; a,; 
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which are trivially related by 



1 ^ 

hi = H + $,, H = -J2h^, (2.8) 

where the dot refers to a derivative with respect to physical time. We define the shear as 

= ^7., (2.9) 

and introduce the scalar shear by 

a^^a,,a'^ =J2P' ■ (2-10) 



(See appendix A of Ref. [23|] to see the relation with the shear usually defined in the 1 + 3 
formalism.) To finish, we define the conformal time by dt = Sdrj, in terms of which the 
metric (12.31) is recast as 

ds^ = S\7]) [-dr/2 + 7,j(r/)dxMx^] . (2.11) 

We define the comoving Hubble parameter by = S' / S, where a prime refers to a derivative 
with respect to the conformal time. The shear tensor, now defined as 

'^ij = ^7;- , (2.12) 

is clearly related to by cTjj = Scij so that cr^ = (Tija^^ is explicitly given by = X]i=i(A')^ 
and is related to its cosmic time analogous by a = Sa. 



B. Friedmann equations and their general solutions 

1. Friedmann equations 

In cosmic time, assuming a general fiuid as matter source with stress-energy tensor 

Tf,u = pUf,Uy + P{g^,^ + n^Uj,) + vr^^, (2.13) 

where p is the energy density, P the isotropic pressure and tt^j, the anisotropic stress {n^^u^ = 
and TT^ = 0), the Einstein equations take the form 

3H^ = ^p+^a\ (2.14) 

I = -^(p + 3P)-ia^ (2.15) 

(aj)- = -3Ha'j + ^r., (2.16) 
and the conservation equation for the matter reads 

p + 3H{p + P) + a.jr^ = , (2.17) 
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where the ij-component of vr^i, has been defined as S^Tiij (so that itj = Y^kj)- 

In this work, we focus on infiation and assume that the matter content of the universe is 
described by a single scalar field so that 



This implies that 



3H^ = K 
S 

S ^ 



-a 



[0^ - viv)] - -a' 



0- 



The last of these equations is easily integrated and gives 



^3 



where /C* is a constant tensor, (/C*)' = 0. This implies that 



a 



with /C^ = /C*/C^, from which we deduce that 



^6 



a = -3Ha . 



(2.18) 

(2.19) 

(2.20) 
(2.21) 

(2.22) 

(2.23) 
(2.24) 



The conservation equation reduces to the Klein-Gordon equation, which keeps its Friedman- 
nian form, 

(p + 3if0 + = 0. (2.25) 

2. General solutions 

Let us concentrate on the particular case in which ti^^ = (relevant for scalar fields) and 
first set 

A = B,W{t) , (2.26) 
where Bi are constants yet to be determined. Equations (12.101) and (12.231) then imply that 

Xa?)*"(*' = ^' 



from which we deduce that 



W{t) 



dt 



The constraints (12.51) and (12.101) imply that the Bi must satisfy 



(2.27) 



(2.28) 



i=l 



1=1 
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which can be trivially satisfied by setting 



/ 2 Stt 
-Bj = A /-/C sin Oj, with ai = a + —i, zG {1,2,3}. (2.29) 
V o o 



Thus, the general solution is of the form 



p.(^t) = j'^ICsmia+'^i] xW{t), (2.30) 



3 V 3 
where S is solution of 

3i7^ = Kp + l^. (2.31) 

Once an equation of state is specified, the conservation equation gives p[S] and we can solve 
for S{t). 

As we shall see, it is convenient to introduce the reduced shear 

1 (7 1 0" , , 

in terms of which the Friedmann equation takes the form 

(1 - x')H' = |p , 
so that the local positivity of the energy density implies that < 1. 



3. Particular case of a cosmological constant 

First, let us consider the case of a pure cosmological constant, V = constant and ip = 0. 
This case is relevant for the initial stage of the infiationary period since we expect first the 
shear to dominate and the field energy density to be dominated by its potential energy. If 
not, then the field is fast rolling and its energy density behaves as exactly as the square 
of the shear. 

The Friedmann equation now takes the form 



1 + 



Q \ 6 



s 

with^ Vo = kV/3 and = (/C^/6Vo)^^^. It can be easily integrated to get 

^(t) = [sinh(t/r.)]'/' (2.33) 
where we have introduced the characteristic time 

= . (2.34) 



1 Note the dimensions: p ~ M^, G - M^^ H ^ M, a ^ M, Bi ^ M,W ^ M~\ Vq ~ JVP, - M° 
where M is a mass scale. 



8 



Thus, we obtain from Eq. (12.271) 



Wit) = Wo + ^\og 



tanh I 

2r, 



(2.35) 



where the constant Wq can be set to zero (corresponding to the choice of the origin of time) . 
It follows, using that a/3/2/C = S'^/r*, that the directional scale factors behave as 





1/3 




sinh ( — ) 




tanh 1 1 






. V2T-J. 



I sin Oi 



(2.36) 



From this expression, we deduce that the directional Hubble parameters evolve as 

1 1 



hi 



[2 sinaj + cosh(t/r^, 



3r* sinh(t/r^, 

and further that the average Hubble parameter and reduced shear are given by 

1 1 1 



H 



3t^ tanh(t/r=^) ' 



X 



cosh(t/r*) 



(2.37) 



(2.38) 



The behaviours of the directional scale factors are depicted on Fig. [T] for various values 
of the parameter a. We restrict to a G [0, 27r/3] and it is clear that for a < n /2 the i = 3 
direction is bouncing. At a = 7r/2 none of the direction is contracting and then it switches to 
direction i = 2, when a > 7r/2. It is thus clear that there is always one bouncing direction, 
except in the particular case in which a = tt/2. 

As a first conclusion, let us compute the time at which the contracting direction bounces. 
Fig [2] (left) depicts the value of the time of the bounce as a function of a and we conclude 
that it is always smaller than IAt^. 



4- Behaviour close to the singularity 



Whatever the potential chosen for the inflaton, the Friedmann equation will be dominated 
by the shear close to the singularity, so that we can use the solutions obtained in the case 
of a pure cosmological constant for the sake of the discussion. From the previous analysis, 
we obtain that 



X; 



t 

2^ 



1 H (1 — sinttj 

18^ 



O 



(2.39) 



The metric can thus be expanded around a Kasner solution of the form 

3 / , N 2pi 



dst 



1=1 ^ ' 



(2.40) 



with the indices 



Pi[a) = -smoj + - 



(2.41) 
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-2.0 -1,5 -1.0 -0.5 0.0 0.5 



FIG. 1: Evolution of the scale factors according to the value of the parameter a. We depict the 
three directional scale factors and the average scale factor S (dashed line) , all in units of 5* . The 
three directional scale factors are permuted when a is changed by 27r/3. Note that there exist 
two particular cases in which the spacetime has an extra rotational symmetry when a = 7r/6 or 
a = 7r/2. The latter case is even more peculiar since this is the only Bianchi / universe for which 
none of the direction is bouncing. 



that clearly satisfy J2Pi = JZpI = ^ [see Fig. [2] (right)]. We thus have 



1=1 ^ ' 



l + -(l-p,) - 



12 



(dx' 



i\1 



(2.42) 



up to terms of order (t/T*) 
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Kasner exponents 





FIG. 2: (left) Time at which the direction i bounces (blue: i = 3), (red: i = 2). At 7r/2, none of 
the direction is contracting and a change of the contracting direction occurs, (right) Value of the 
Kasner exponents as a function of the parameter a characterising the Bianchi / model close to the 
singularity. 



The invariants of the metric behave as 

4 

3^ ' 
4 

9^ ' 
1 



R 



RfiuR^'^ 



+ 



4 



cosh' \t/{2n)] 
32 cosh {t/% 



sinh [t/r^ 



[3 



+ 



COS a sm a 



sm a 



1] 



(2.43) 
(2.44) 



(2.45) 



Clearly, we see that R and R^^R^'^ are regular at the singularity, which is expected for a cos- 
mological constant dominated universe since R^y = —h^g^j^. The third invariant Rf^i^p^R'^'^'"^ 
diverges when we are approaching the singularity, the only exception being the case where 
a = 7r/2 which corresponds to the positive branch considered in Ref. 2J]. In this particular 
case, the Kasner metric has exponents (1,0,0). We emphasize that the spacetime a = 7r/2 
is a singular point in the set of Bianchi spacetimes since there is no uniform convergence of 
the invariants of the metric evaluated on the singularity when a ^ 7i/2. 



C. Slow roll parameters 

In order to discuss our results, we introduce the slow-roll parameters in the usual way by 

n'2 



6 



ip''^ + 2S^V' " ' Hip' Hif 
With these definitions, the Friedmann and Klein-Gordon equations take the form 



;i - x')n^ 



K 



and we deduce that 



3-e 
W 



-VS' 



(3 - 5)nip' + V^S^ = , 



;i-6) + (e-3)x^ 



(2.46) 

(2.47) 
(2.48) 
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so that 

^ = n^[2-e + {e-3)x''] (2.49) 

and 

e = 2ne{e - S) . (2.50) 

Interestingly, it is easy to check that 

x' = — 7Yx(l — x^)(3 — e) , 

from which we deduce that 



6' = n 



-9x' + - (1 - x^)i3e + 35) + 6(6 + e(l - x^)) 

rC 



(2.51) 



From the definition f l2.46p and making use of the second equation of Eqs. fl2.47p . we 
deduce that the slow-roll e parameter takes the form 



2k \V J \?,-6 
Once we use the relation tp^ = 2eV/ (3 — e), we deduce that 



As long as 0^ <^ ^i'^) and in the limit a; — >• 1, which corresponds to the shear dominated 
period prior to the infiationary phase, we have that 

e ^ , 5^-3, (2.53) 

but we still have 6' 0. It follows that initially, even if the shear decreases rapidly, ip 
remains almost constant and 6 remains close to —3. This solution converges when t — to 
the pure cosmological constant solution of §IIB 31 Note that this conclusion differs from the 



statement of Ref. 39 



Before, we reach the slow-roll attractor, we may be in a transitory regime in which 
0^ ^ V. Then, this implies that 

e 3, 5 -> 3, V-Y' 

The field velocity decreases so that this solution converges rapidly toward the slow-rolling 
attractor. These general results on the limiting behaviours will be important to understand 
the dynamics of the infiaton. 



D. Numerical integration for a massive scalar field 

In this article, we will consider the explicit example of chaotic infiation with a potential 

V = JmV^ • (2.54) 
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The dynamical equations (12.19p and fl2.25p can be rescaled as 



6 



(2.55) 

ip + 3hip + ip = , (2.56) 



where we use r = mt as a time variable (so that h = H/m), ip = (p/Mp (with ^ = 
8ttG = k) as the field variable, and where 5** = (/C/m)^/"^. Under this form, it is clear that 
the various solutions of the system fl2.55H2.5^ are, in general, characterised by the three 
numbers {ipito), ip(to), S^:}. As we will now see, this extra dependence on the parameter 5"^, 
causes our dynamical system to behave differently from its analogous in Friedmann-Lemaitre 
spacetimes^. Under this form, we also see clearly that one cannot continuously go from a 
Bianchi to a FL spacetime (because either S** = or 7^ 0), even though Bianchi spacetime 
isotropizes. 

It is well known that the dynamics of the inflationary stage in a Friedmann-Lemaitre 
spacetime is characterised by attractor solutions, clearly seen in the phase space. For a large 
set of initial conditions, the solutions converge to the slow-roll stage deflned by an almost 
constant (f and (f decreasing accordingly. It is thus given by a roughly horizontal line in 
the phase portrait which ends with the oscillations of the scalar fleld at the bottom of its 
potential. As previously mentioned, the slow-roll inflationary stage in Bianchi I spacetimes 
possesses the same attractor behaviour, although it is quite different during the initial shear- 
dominated phase. When the shear dominates, the solutions are rapidly attracted to the point 
~ and if nearly constant. This is the attractor whose solution is given by Eqs. (12.581 
I2.62p . It then converges towards the Friedmann-Lemaitre behaviour when the shear becomes 
negligible (see Fig. H]). 

In conclusion, we have a double attraction mechanism, namely of the field dynamics 
toward the slow-roll attractor and of the Bianchi spacetime toward a Friedmann-Lemaitre 
solution. 

In Fig. \5[ we compare the dynamics of the infiationary phase of a Bianchi / and a 
Friedmann-Lemaitre spacetimes. While in both cases the (average) scale factors grow, the 
effect of a shear-dominated phase is to initially decrease the velocity of the expansion of the 
Bianchi universe. On the other hand, the dynamics of the scalar field is barely affected by 
the presence of the shear as long as we have reached the slow-roll attractor. In the following 
we will assume that we have reached this attractor. If the number of e-folds is small, it is 
not clear that this attractor has been reached, regardless of whether or not we assume a 
Friedmann-Lemaitre spacetime. In such a case, one has to rely on families of trajectories to 
draw the observational predictions (see e.g. Ref. 4o|). 



Let us now consider the behaviour close to the singularity and introduce the characteristic 
time r*, explicitly given by 

2" M„ , , 

. 2.57 

6 rrupQ 



^ There is another way to see this: the term in in Eq. (|2.55p is equivalent to the contribution of a 
massless scalar field, x say: that would satisfy the Klein-Gordon equation x + S/ix = 0, implying x S^^. 
One requires initial conditions for this extra-degree of freedom as well. 
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versus if (FL) 



If versus ip (BiaQchi) 





FIG. 3: Comparison of the phase portraits of a Friedmann-Lemaitre (left) and a Bianchi / infla- 
tionary phases. The curves with same colour corresponds to same initial conditions for the scalar 
field but starting with different initial shear (see Fig. [J] for a 3-dimensional representation) . 




FIG. 4: Phase portrait of a Bianchi / inflationary phases in the space {(p,ip,a}. The plane a = 
corresponds to the FL-limit (see Fig. [3]). This illustrates the double attraction mechanism of the 
spacetime toward FL and of the solution toward the slow-roll attractor. 
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Log[mt] 



Log[mt] 



FIG. 5: Evolution of the average scale factor (left) and scalar field (right) as a function of time 
for both the Bianchi (red, solid line) and FL cases (blue, dashed line). The Bianchi solution is 
characterised by ipo = 16Mp. We have normalized the solutions such that the scale factors have 
the same value at the end of the inflation when the shear is negligible. The inner left figure shows 
the logarithmic evolution of the velocity of the expansion in both cases, the minimum of which 
indicates the time at which S = for the Bianchi case. 



By developing all the previous equations in powers of t/r^,, we obtain that 



Hit) 
xit) 

m 

e{t) 



1 

3t 



1^ 

3r2 



If 



1-7:^ + — 



2t: 



6 Vv^o 

2^ + 



o 



O 



3 V '/'o 



(2.58) 
(2.59) 

(2.60) 
(2.61) 
(2.62) 



in complete agreement with the expansions obtained in Ref. [2^. In this limit, we understand 
why 5 —3 close to the singularity. It simply reflects the fact that the field decreases as 
. This is different from the case of chaotic inflation in a Friedmann-Lemaitre spacetime 
where the field varies linearly with time during the slow-roll regime since in that case 



^_fMX t_ 

\Vi } T* 



S = Si exp<j jpi^'^i -v'^it)] 



and the slow-roll parameters are explicitly given by 

M2 



2^ 



6 = 



Let us come back on the slow-roll parameters defined in § III CI In the particular case of 
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a quadratic potential, we have 



5' = n{3 - 5) 



36 



x\3 



(2.63) 



since S'^V^^ = eTi}{3 — 5)^/(3 — e). During the shear-dominated phase, 5 ~ — 3 and the 
previous equation tells us that 5 remains constant. It follows from Eq. (12.501) that e' is of the 
same order than e (note the difference between the standard case in which e' is second order). 
It follows that the variation of e cannot be neglected until 5 has converged toward 0. Then e 
can be considered as constant until the end of inflation. While the universe isotropizes, both 
5 and e converged toward their Friedmann-Lemaitre value. Fig. [6] illustrates the evolution 
of the two slow-roll parameters and compare them to their values in a Friedmann-Lemaitre 
universe. 





FIG. 6: Evolution of the slow-roll parameters e (left) and 5 (right) during the inflationary phase 
for a Bianchi / model (solid lines) with a = tt/A compared to the case of a Priedmanian model 
(dashed lines). We assume a potential of the form (j2.54p and initial values ipo = (16,26,36)Mp 
corresponding respectively to the red, green and blue lines. 

For any initial value of the scalar field, ipo, we define ipi as the value of the field when the 
universe starts to inflate, that is at the time when 5" = 0. Then, the number of e-folds of 
the inflationary period is defined as 

A'M.bf^^^ r'%. (2.64) 

where is the value of the field at which e = 1. Since at that time the shear is negligible, 
is given by 

= V2Mp . (2.65) 

It follows that for any initial value of the field, we can characterise the shear-dominated 
phase by 

Aipl^po] = ifo- v^ibo] , (2.66) 
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which indicates by how much the field has moved prior to infiation. Then, the duration of 
the infiationary phase is given 



(2.67) 



This has to be compared with the number of e-folds in the FL case 

2 

I / ID: \ 



Mr, 



1 

2 



In Fig. [7] (left), we depict the fractional duration of the shear-dominated phase [see Eq. 
fl2.66p ] as a function of (fQ. We see that the larger the initial value of the field, the smaller 
the impact of the shear. Also shown in Fig. [7] is the number of e-folds for the Friedmann- 
Lemaitre and Bianchi spacetimes. The presence of the shear slightly increases the number 
of e-folds (solid red fine in Fig. [7]). 
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FIG. 7: (left) Relative variation of the scalar field during the shear dominated era in function of 
its initial value, (right) Comparison of the number of e-folds for a Bianchi (solid, red line) and 
Friedmann-Lemaitre (blue, dashed line) inflationary period as a function of ^pi. The inside plot 
is the relative difference between the number of e-folds of the Bianchi and Friedmann-Lemaitre 
inflationary period as a function of <pi. 



E. Discussion 



This general study of the background shows that generically there is always one bouncing 
direction, except in the particular case a = 7r/2 considered in Ref. [1^ (positive branch). 
Regarding the dynamics of the universe, the expansions (12.58112.6^ at lowest order exactly 
reproduce the exact (numerical) solutions. 

This analysis actually shows that the shear is effective typically until the characteristic 
time r*. Going backward in time, the universe goes rapidly to an initial singularity and is 
thus past incomplete (in fact as other inflationary models (4l|)- Indeed, as we shall discuss 
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later, we do not want to extrapolate such a model up to the Planck or string times and we 
just assume that they are a good description of the inflating universe after this time. 

We have focused on the evolution of the "slow-roll" parameters and concluded that the 
variation of e cannot be neglected until the shear has decayed. We have also shown that the 
scalar field barely moves during the shear dominated era and that, with same initial value 
for the scalar field, the number of e-folds is almost not affected (even though slightly larger) 
by a non- vanishing shear. 



III. SUMMARY OF THE PERTURBATION THEORY 

In our previous work , we investigated the theory of cosmological perturbations around 
a Bianchi I universe. We shall now briefly summarize the main steps which are necessary 
to deduce the dynamical system of equations that we want to solve in this article. 



A. Mode decomposition 

First, we pick up a co moving coordinates system, {x*}, on the constant time hypersur- 
faces. Any scalar function can then be decomposed in Fourier modes as 

/(^^^)= / ^^2nk.v)e'''''. (3.1) 

In the Fourier space, the comoving wave co-vectors fcj are constant, k'- = 0. We now define 
/c' = that is obviously a time-dependent quantity. Contrary to the standard Friedmann- 
Lemaitre case, we must be careful not to trivially identify ki and /c*, since this does not 
commute with the time evolution. Note however that Xik^ = x^ki remains constant. 

Besides, since (/c*)' = —2a'^^kp, the modulus of the comoving wave vector, k'^ = k'^ki = 
Y-'kikj, is now time-dependent and its rate of change is explicitly given by 

J = -a'^kkj , (3.2) 

where we have defined the unit vector 

k^^. (3.3) 

k 

Now, we introduce the base {e^,e^} of the subspace perpendicular to k\ By construction, 
it satisfies the orthonormalisation conditions 

e^k,f= = , e^e'^^'^ = 5^' . 

Such a basis is indeed defined up to a rotation about the axis k\ These two basis vectors 
allow to define a projection operator onto the subspace perpendicular to /c* as 

Pij = e\e) + e\e] = %j - hk ■ (3-4) 

->i pi 



It trivially satisfies P'^Pl = P^k^ = and P'^-fij = 2. 
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Any symmetric tensor Vij that is transverse and trace-free has only two independent 
components and can be decomposed as 

V^,ih,v)= E Vxik\r^)ef^{h) , (3.5) 

A=+,x 

where the polarization tensors have been defined as 

p'^ p'^ p^ p'^ _|_ p2 p 1 

4 = "'"^J'"^ 6l + ^^^^^5^ . (3.6) 

It can be checked that they are traceless {efjY'' = 0); transverse {efjk^ = 0), and that the 
two polarizations are perpendicular {Sij^^/i = 5^)- This defines the two tensor degrees of 
freedom. 



B. Decomposition of the shear 

It is then fruitful to decompose the shear in a local basis adapted to the mode that we 
are considering. The shear being a symmetric trace-free tensor, it can be decomposed on 
the basis {ki, ej, e|} as 

CTrj = ^ (kkj - ^7ii ] ^11 + 2 X] f^va %e") + E (T^xSij- (3.7) 
^ ^ a=l,2 A=+,x 

This decomposition involves 5 independent components in a basis adapted to the wavenum- 
ber ki. We must however stress that (cTh , (Xya, o"-pa) do not have to be interpreted as the 
Fourier components of the shear, even if they explicitly depend on ki. This dependence 
arises from the local anisotropy of space. 

Using Eq. (13. 7p . it is easily worked out that (Tjj7*-' = 0, and that 

aijk = a^^kj + E cTvae" , aijky = a^^ , (3.8) 

a 

and 

aijs'l = a^x , (^ijk'ei = a^a ■ (3.9) 
The scalar shear is explicitly given by 

3 



a.,a- = V + 2j]a;, + j:aJ,, (3.10) 



A 



which is, by construction, independent of ki. We emphasize that the local positivity of the 
energy density of matter implies [see Eq. ( 12.32^ ] that < H? and thus 

-a„ < -^a <n. (3.11) 

2 II - V6 

This, in turn, implies that 

< 2n, (3.12) 

a property that shall turn to be very useful in the following of our discussion. Analogously, 
we have that 

a^A<V6^- (3.13) 
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C. Gauge invariant variables 

We start from the most general metric of an almost Bianchi / spacetime, 
ds^ = S^[-{1 + 2 A) + 2Bidx'dr] + (7^^- + hij) dx'dx^] . 
Bi and hij are further decomposed as 

Bi = diB + Bi , 



with 



diB' = = diE\ Ei = = diE'^. 
We showed, that one can construct the following gauge invariant variables 

{k^E)' 



^ = A+^\s B 



^) = -C~H B- 



jk'E)'^ 
k^ 



for the scalar modes, 



= - 7i,- {E^y + 2ik^aijP[E, 



(3.14) 

(3.15) 
(3.16) 

(3.17) 

(3.18) 
(3.19) 

(3.20) 



for the vector modes and that the tensor mode Eij is readily gauge invariant. 

Concerning the matter sector, one can introduce a single gauge invariant variable associ- 
ated with the scalar field perturbation. 



Q = 6ip- ^if' . 



(3.21) 



D. The Mukhanov- Sasaki variables and their evolution equations 

We established that the only degrees of freedom reduce to a scalar mode and two tensor 
modes 

v = SQ, v^/iA = SEx . (3.22) 

They evolve according to 



v" + ujl{ki, r])v = ^ ^x{ki, r])nx , 

A 

+ ujl{ki, 'r])fix = ^x{ki, v)v + ^(fci, ?7)/i{i-A) 
where the pulsations are explicitly given by 



2 . .2/1^ „\ — l2 



(3.23) 
(3.24) 

(3.25) 
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The two functions z^^ and z\ have been defined by 



/ 



+ . (3.26) 



^ ? + 2aJ„_„ + 4fsV,)' + 4(l^l . (3.27) 



/ 



and we also need the couphng terms 



1 ^ / 



Mv,h) = i^v^ 1^ 7^17' 1 ' (3-28) 
^iv, h) = ;^ ( ^^I^ 1 - 2^TX^.+ . (3.29) 



IV. PRESCRIPTION FOR THE INITIAL CONDITIONS 

The set of equations (13 .23113.291) completely determines the evolution of the three degrees 
of freedom of our problem. To be predictive, we must determine the initial conditions. 

In a FL spacetime, the procedure is well understood and relies on the quantization 
of the canonical variables on sub-Hubble scales where it can be shown that they evolve 
adiabatically. 

We have to understand how far this procedure can be extended to a Bianchi universe 
and how robust it is to the existence of a non- vanishing primordial shear. We thus start, in 
§ IIV Al by a review of the standard FL procedure, to highlight its hypothesis. In § lIVBt we 
stress, and also quantify, the differences that appear in a Bianchi universe. This will lead 
us (§ lIVCp to propose an extension of the quantization procedure. We shall finish in § IIVDI 
by critically discussing the limits and weaknesses of our quantization procedure. 

A. Friedmann-Lemaitre universes 

For simplicity, let us consider the case of a pure de Sitter phase. This represents no 
limitation to our following arguments and can be generalised to an almost-de Sitter phase. 

1. Quantization procedure 

In order to be quantized, the canonical variables are promoted to the status of quantum 
operators jl] and are decomposed as 

f d^k r 

= j [t^fc(r?)e'''-"ak + t;^(r?)e-*'^-X 

^^[i)k(r/)e^'^- + 4(r^)e-''-] , (4.1) 
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where the creation and annihilation operators satisfy the commutation relations [cik, aj^/] = 
5'^^)(k — k'). The mode function, Vk{f]), is solution of the classical Klein-Gordon equation 

vl + iul{k,ri)vk = with u;l{k,ri) = e-^ , (4.2) 

which is the equation of motion for a harmonic oscillator with time-dependent mass, which 
translates the fact that the field lives in a time-dependent background spacetime. The 
general solution of Eq. (14.21) is 

Vkiv) = [Aik)Hi^\-kr]) + B{k)Hi^\-kr])] ^ 

where Hy are the Hankel functions. In the particular case of a de Sitter era considered here, 
= 3/2 so that 




TTZ \ IZ 



Canonical quantization consists in imposing the commutation rules ['0(x, 77), ^(x', 77)] = 
[7r(x, 77), :n"(x', 77)] = and [?}(x, 7/), 7r(x', 7/)] = 5*^^''(x — x') on constant time hypersurfaces, 
Tt being the conjugate momentum of v. From Eq. (14.11) and the commutation rules of the 
annihilation and creation operators, this implies that 

Vkv'^ - vlv't, = i , (4.3) 

which determines the normalisation of the Wronskian. The choice of a specific mode function 
'^kiji) corresponds to the choice of a prescription for the physical vacuum |0), defined by 

SklO) = . 

The most natural choice for the vacuum is to pick up the solution that corresponds adia- 
batically to the usual Minkowsky vacuum so that 

Vk ^ ^e-^'^'' 



'2k 

when kr] —00. This implies that the mode function is 



'2k \ ikr] 

This choice is referred to as the Bunch-Davies vacuum 



v, = —={l + -L) e-^'^^ . (4.4) 



2. WKB approximation 

In more general cases, and for sure in the Bianchi case that follows, we may not have 
exact solutions for the mode functions. On can redo the previous construction by relying 
on a WKB approach [12] in which one introduces the WKB mode function 

.r"(7,) = . (4.5) 
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It is easily checked that it is solution of 

„,WKB" j_ (, ,2 ^ \ „,WKB _ r, 

with 

QWKB = T — - o— • (4-6) 

For a function satisfying an equation such as Eq. (14. 2p . the WKB solution is thus a good 
approximation as long as the WKB condition IQwkb/i-^^I -C 1 is satisfied. In the example at 
hand, this condition reduces to /c?7 — — oo so that on sub- Hubble scales, the mode function 
is actually close to its WKB approximation. We see that the quantization procedure thus 
relies on the fact that there exists an adiabatic (WKB) solution on sub-Hubble scales. 



3. Primordial spectra on super-Huhhle scales 

Once the initial conditions are fixed, Vk is completely determined and it can then be 
related to the scalar field perturbation Q (also promoted to the status of operator). After 
the modes become super-Hubble, i.e. kr] <^ 1, the scalar field perturbation in flat slicing 
gauge is given by 

where we have used that S{ri) = —l/Ht] for a pure de Sitter inflationary phase, H being a 
constant in this case. All the modes are proportional to (ok + al_k) ^'^ that the variables Qk 
commute. We thus deduce that Q has actually the same statistical properties as a Gaussian 
classical stochastic field. Effectively, we can replace our quantum operators by stochastic 
fields with Gaussian statistics and we introduce a unit Gaussian random variable, et,(k), 
which satisfies 

(e.(k))=0, (e.(k)e:(k')) = 5(^)(k-k') . 

In this description the mode operators are replaced by stochastic variables according to 
^'k ^ f k = Vk{r])ey{'k) and we identify the (quantum) average in the vacuum, i.e. (0|...|0) by 
an ensemble (classical) average, (...). 

The correlation function of v is defined as 

= (0|i)(x,r/)i)(x',r/)|0) , 

and takes the simple form 

Interestingly, in a Friedmann universe, isotropy implies that we can integrate over the angle 
to get 

f dk , ,n sin kr , , 

and it is , because of the symmetries of the background, a function of r = |x — x'| only. We 
thus define the power spectra 

Pvik) = \vk\'' , r,{k) = ^\vk\^ . (4.9) 
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In the stochastic picture, the correlator of is simply given by 

{vkvl,) = P,{k)6'^'\k-k') , (4.10) 
from which one easily deduces the power spectrum of the curvature perturbation 



Indeed, we can proceed in the same way for gravity waves. Since they are not coupled 
to scalar modes and since the two polarisations are independent, we introduce two sets of 
creation and annihilation of operators, 6k, a, one per polarisation. On super- Hubble scales, 
the two modes can be described by two independent Gaussian classical stochastic fields, 
/ik,A = /ikeA(k) with 

(eA(k)el,(k'))=(5AA'5(^nk-k'). 
The power spectra are thus given by 

27r2 

Px{k) = -^^xik) = l/ifcp . 
Spatial isotropy implies that P+ = Px so that the tensor modes power spectrum is 

'S 



Prik) = 2-P^{k) . (4.11) 



^. Conclusion 



To conclude, this short review of the standard procedure highlights (1) the importance 
of the WKB regime on sub-Hubble scales which allows to construct a Bunch-Davies vacuum 
adiabatically, (2) the fact that on super-Hubble scales (where one wants to draw the predic- 
tions for the initial power spectra) the quantum operators can be conveniently replaced by 
stochastic fields, (3) the importance of isotropy which implies that there exist 3 independent 
stochastic directions (because modes are decoupled) and (4) the fact that the two gravity 
wave polarisations have the same power spectrum. 

We shall now see which of these properties generalise to a Bianchi universe. 



B. Generic Bianchi I universes 



1. Characteristic wavenumber 

In order to relate our predictions to observations, we introduce the characteristic 
wavenumber /cref by 

kref ^ SH\,^^^ . (4.12) 

We define M as the number of e-fold with the quantity SH rather than S in the defini- 
tion 02.641) . If we denote by Afrd the number of e- folds between t = and the end of 
inflation, and A/q the number of e-folds from the end of inflation until now, then we can 
relate k^^i to the largest observable scale today, fco = SqHq, by 

^ = gCATo-Met) . (4.13) 
ko 
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Since during the inflationary era after t^,, H is nearly constant, then ACrf — A^ref- As for TVo, 
it depends on the post-inflationary evolution and it can be estimated 4^ by 



. 62 - ,„ ( ™) . i 1. ^ - i 1. (§) - .n K (4.14, 

h being the Hubble parameter in units of 100 km.s~^/Mpc. Typically, Vk^ ~ Vend as long 
as slow rolling holds. The reheating temperature can be argued to be larger than pl^^ > 
10^° GeV to avoid the gravitino problem [4J]. The amplitude of the cosmological fluctuations 
(typically of order 2 x 10~^ on Hubble scales) roughly implies that V^n(f is smaller than a 
few times 10^^ GeV and, for the same reason as above, has to be larger than 10^*^ GeV in 
the extreme case. This implies that Nq has approximately to lie between 50 and 70, which 
is the order of magnitude also required to solve the horizon and flatness problem. 

iVref can be computed from the background dynamics, and as shown on Fig. [3, it is almost 
equivalent to its value in the FL case. 



2. Anisotropy 

The first obvious difference with the FL case arises from the local spatial anisotropy. 

First, it is clear from the set of equations (13.23113.2^ that it implies that the evolution 
of the mode functions shall depend on ki and not simply on the modulus. This violation of 
isotropy will reflect itself on the fact that 

1. the power spectra at the end of inflation will be functions of k and not k, i.e. Pv{ki), 

2. because of the coupling between scalar and gravity waves, there exists a cross- 
correlation between scalar and tensor, i.e. {vfix) ^ 0; 

3. the two polarisations shall a priori have two different power spectra, i.e. P+ 7^ P^. 

A second related issue arises from the evolution of a comoving wavenumber. Let us con- 
sider the different evolutions of a wave-mode of modulus k at the end of inflation according 
to its orientation. As we see on Fig. [SI depending on its orientation, this mode has very 
different time evolutions before it settles to a constant value. 



3. WKB regime 

Whatever the Bianchi universe we consider, there is a shear-dominated phase prior to 
inflation. During this phase S < 0. Besides, in a generic Bianchi I spacetime (that is 
a ^ vr/2) two of the scale factors go to zero while the third is bouncing (in the a = 7r/2 
case, one scale factor goes to zero while the two others remain constant). 

For analysing the WKB regime we will consider, as usual, the ratio k/SH to discuss 
whether a given mode is inside {k/SH ^ 1) or outside {k/SH ^ 1) the Hubble radius. 
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FIG. 8: Logarithm of the ratio k/k^-ef where /cref is the modulus of the wavenumber when the shear 
is zero. We see that as soon as the shear grows, k depends on the direction on which it is ahgned 
(Each color corresponds to a principal axis of the Bianchi universe). We have considered a generic 
Bianchi universe with a = 7r/4. 

According to the details shown in the Appendix A, if we specify to only one direction we 
have /c ~ l/cj ~ S/Xi, then 

^ = ^ ^ .2{l-sinai)/3 

SH XiH 

during the shear-dominated regime. This shows that, except when a = 7r/2, any given mode 
becomes super- Hubble when we approach the singularity (t — > 0), and that this approach is 
faster (going backwards in time) for modes aligned with the bouncing direction. 

In other words, all modes become super-Hubble in the past, and the mode aligned with 
the bouncing direction (blue line in Fig. [S]) becomes super-Hubble earlier (again, going 
backwards in time) than the ones with same k at the end of inflation but aligned with a 
growing direction (blue line in Fig. [H]). 

For these two reasons, we can doubt the existence of a well-defined adiabatic vacuum 
for all modes through their early evolution, as happens in FL universes. However, we can 
still ask whether the WKB regime is reached in a short time interval when the shear is not 
complete negligible, and how long it lasts given the wavenumber. 

Let us thus discuss quantitatively the validity of the WKB approximation. First, we focus 
on the pulsation, we neglect the effect of the couplings, and consider the WKB solutions 

.r-(,) = 1 e^^/-^^ , ^^Z''iv) = -^e±^/-<^^ . (4.15) 

They are good approximations of the solutions of Eqs. (13. 23113. 24|) if IQ^x'^ / uj"^ x\ ^ ^ where 
QwKB has been defined in Eq. (14. 6p . 

Fig. [9] illustrates the validity of the WKB approximation for the scalar modes. It depicts 
the evolution of IQ^^^^/u^l as a function of time for three different modes corresponding 
to the three principal axis of the Bianchi universe. We see that, for a given comoving 
wavenumber k at the end of inflation, the WKB condition is always violated in the past and 
that it is violated first in increasing order of the Kasner coefficients (compare with Fig. [8]). It 
can be checked that the larger the k is, the longer the time during which the WKB condition 
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is restored. For long wavelength modes (typically of order l//crcf), the WKB regime is never 
established. 

Fig. [To] compares the exact (numerical) solution and the WKB approximation for a mode 
k = lOfcref- Fig. [TT]is similar to Fig. Obut for the tensor modes. Indeed, we reach the same 
conclusions. 
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FIG. 9: Left panel: Evolution of the quantity |Q/a;^| for a FL universe. Right panel: Evolution 
of IQ^^^/w^l for three different modes, each of them aligned with one of the three orthogonal 
directions (same color code as in Fig. [8|), and with the same modulus lO/crcf at the end of inflation. 
The vertical line correspond to the instant r/(T*) and we have considered a generic Bianchi spacetime 
with Q = 7r/4. 




FIG. 10: Comparison of the exact (solid lines) solutions and the WKB (dashed lines) solu- 
tions (j4.15p for the scalar modes for three different directions with k = lOfcref at the end of 
inflation. We have considered a generic Bianchi spacetime with a = 7r/4. The figures show two 
genuine WKB modes (left panel) and one non-WKB mode (right panel). The inner-left figure 
shows the ratio |t^^^^/uk| for the modes which satisfy the WKB approximation. 



4- Couplings 

The third difference arises from the coupling between scalar and tensor modes. 
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FIG. 11: Evolution of IQ^^^/w^l for the two tensor polarisations (left: A = +, right: A = x) and 
for various modes with the same modulus k = lOfcrcf at the end of inflation. The vertical lines 
represent the time i](t^,). We have considered a generic Bianchi spacetime with a = 7r/4. 



As we demonstrated in our previous analysis j23[ (see § IV.D), deep in the sub-Hubble 
regime (that is when k/SH is large enough), the three degrees of freedom decouple and 
behave as a collection of three independent harmonic oscillators. 

However, on larger scales the couplings are a priori non-negligible. We thus need to 
evaluate with care the scales for which this is a good approximation. 

Fig. [T2] shows that, while the functions Ka, which couple gravity waves and scalar modes 
can be neglected on small scales, this is certainly not the case for the coupling H between 
the two gravity wave polarizations. This coupling cannot be neglected, even at early time, 
for modes which are not sub-Hubble enough. Typically, the modes for which this coupling 
cannot be neglected correspond to modes for which the WKB regime cannot be reached. 

Let us compare this situation with the case of multi-field inflation. The equations of 
evolution for the various scalar field perturbations are also usually coupled (see e.g. Ref. 45 
for a recent review). However generally, it is possible to extract independent fields, at least 



in the sub-Hubble regime j46| so that one can introduce a set of independent stochastic 
fields. To our knowledge, the situation where this is not possible has not been addressed. 

The situation is analogous for us and the long wavelength modes will be particularly 
difficult to trustfully deal with because they can never be considered as independent. 




FIG. 12: Evolution of ^x/H? (left: A = -h, middle: A = x) and l/Ti'^ for three modes each of 
which is aligned with one of three orthogonal arbitrary directions (represented by three different 
colors). We have considered a generic Bianchi spacetime with a = 7r/4. Note that these functions 
depend only on the direction of the wavenumber and not on its modulus. 
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C. Prescription for setting the initial conditions 



1. Prescription for Bianchi I spacetimes 

It is part of the nature of quantum fluctuations that they do not have initial conditions 
in the sense that they are continuously excited. As soon as a mode can oscillate, it will 
be sourced by these quantum fluctuations. Given the discussion of the preceding section, 
we will thus set the initial conditions at the time where a mode is the deepest in the WKB 
regime. This time, ti(k) say, depends explicitly on the mode, which is not a problem since all 
the modes are independent from each other. In practice, we have fixed ti(k) by minimizing 

As we saw, modes with k < k^ef never enter a WKB regime. Given the fact that modes up 
to approximatively /cq have been excited and have, at least as a good approximation, a scale 
invariant power spectrum, we have to assume that k^^i ^ ^o- Then, from Figs. [T^ and [T^ 
we deduce that the WKB is reached for all modes with k > k^ef- Indeed this amounts to a 
fine tuning on the shear such that only the largest observable modes today were affected by 
the Bianchi phase. 

To check the robustness of this procedure, we have varied the time ti(k). In particular, 
we have also assumed, as a test, that ti{k) = for all modes. It can be shown that this 
does not affect the predictions for the modes with k > 2/crcf while long wavelength modes 
are more affected. This is simply due to the fact that the duration of their WKB phase is 
smaller (see Fig. [12]). Also note that the procedure is more robust for the two directions 
which are not bouncing. 

For these modes, and as can be seen from Fig. [T51 it is also a good approximation to 
neglect the coupling terms in Eqs. (13.23113.2^ . We emphasize that this hypothesis breaks 
down approximatively at the same time when the WKB approximation also breaks down. 

Therefore we will assume that the three modes are independent at the time when they 
are excited by the quantum fluctuations. 

Technically, we thus start our computation by setting 

'"^ = n N ^^('^) ' ^^'k = , , eA(k) , (4.16) 

up to an arbitrary relative phase which can be absorbed in the definition of the unit random 
variables and where the three random variables satisfy 

{ex{k)e*y{^'))=6xYS^'\k-k') . 



2. Freedom in a time redefinition 



As mentioned in Ref. 24 



the dynamical system of equation for the Mukhanov- Sasaki 
variables admits time reparameterisations that conserve the canonicity of the system, in the 
sense that if v{ri) is a canonical variable that satisfies 



and 



. dv 



drj dr] 
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FIG. 13: Validity of the WKB approximation at the time we set the initial conditions, (left) We 
set the initial conditions at ij(k) and (right) at r*. The quantity \Q^^^ /uj'^,\ is shown for three 
orthogonal modes and we have considered a generic Bianchi spacetime with a = 7r/4. 
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FIG. 14: Validity of the WKB approximation for tensor modes. We show the quantity 
for A = + (left panel) and A = x (right panel) for three orthogonal modes. We compare setting 
the initial conditions at ti(k) (top) and at r^, (bottom). We have considered a generic Bianchi 
spacetime with a = 7r/4. 



then there is a function / and a time r defined as 

f{r,fdT = d'n (4.17) 

through which we can define a new variable u = fv that satisfies 

d^M 2 ^ 1 dw* ^du 

—r + u.n = and u— u — = i . 

dr2 " dr dr 

If, for example, f is a variable for which the WKB condition does not hold, we might wonder 
whether exists a function / that would lead to a different conclusion concerning the validity 
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FIG. 15: Statistical independence at r*. K/A;^ (left) and li/k'^ (right) as a function of /c at t = r^,. 

of the WKB condition for the variable u. 

In the present case, the transformation Eq. (14.171) would lead to the same equations of 
motion satisfied by the new canonical variables 

V = fv, fi+ = fii+, fix = fl^x , (4.18) 
where the new pulsations are defined according to 

and the coupling functions become 

i^A^^, n^-^. (4.20) 

For the sake of simplicity, let us drop and D in the following discussion. It has been 
shown in Ref. 2^ that, if cij„ and ux satisfy the WKB condition, then / is required to satisfy 
the condition 



1 
7 



d"/ 



< LJ^ , with n = 1 . . . 4 (4.21) 



in order for the new equations to also satisfy the WKB condition. Under such conditions, 
it would lead to the same quantization procedure. Note that the condition for n = 4 is 
required when we take \Q^^'^ /uj'^\ -C 1 for the (correct) WKB condition rather than just 



io' jio^ ^ 1, as assumed in Ref. [24 



In our case the pulsations uj^ and uj\ scale like SH when t ^ 0, or like l/r/ in conformal 
time. Let us choose the integration constant Wq of Eq. (I2.35P such that the initial singularity 
corresponds to 77 = 0. Now, if ~ C/?7 then the associated function Qwkb behaves as 



QwKB 



a;2 4C2 

Let us consider now a time redefinition associated to a function /(?]) = 77^. Then the WKB 
condition involves _ ^ 

QwKB + 2) 



- . (4.22) 

C^2 C2 + A(A + 1) ^ ' 

We conclude that it is possible for the WKB condition to be fulfilled with a new time 
coordinate by choosing A = —1/2, provided \C\ 7^ 1/2. 
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Unfortunately, using the expansion (lB2p and the asymptotic behaviours obtained in 
Eqs. f l2.53p . we deduce that when the leading term of is z"/zs, then u"^ — > —H^, which 
implies that for this pulsation, \C\ = 1/2. 

Thus, it is never possible to construct a time redefinition which would enable the WKB 
condition to be satisfied for u^. The only exception arises when the leading term of u"^ is 
k^. This happens for a = 7r/2 since in this particular case k ~ ~ 77"^ and we recover 



the conclusions reached in Ref. 24 



As a conclusion, though we might naively think that redefining time could lead to equa- 
tions satisfying the WKB conditions for the new canonical variables, it is impossible however 
to build such a change of time coordinates. The choice of the canonical variables is thus 
unique up to the reparameterisation satisfying the conditions fl4.2ip . which would not change 
our predictions. 



D. Discussion 



Let us discuss our procedure to set the initial conditions. 



First for modes smaller than /cj-ef , the WKB regime was never reached. We have no natural 
prescription to determine their amplitude. A solution, that we do not investigate in this 
article, would be to fix them by assuming that they minimize their energy, as proposed in 
Ref. 47| in the study of some trans-Planckian models in which the WKB regime is violated. 

On the other hand, and probably in a more conservative way, one could just assume their 
initial value to be completely random and introduce a free function to describe the initial 
conditions on large scales. Such a function would then need to be measured from e.g. large 
angular scale properties of the CMB or predicted by some processes that arise at the Planck 
or string scale, and that indeed cannot be accounted for in our description. 

In the former case, we loose the predictive power on large scales, that actually may just 
be beyond the actual size of the observable universe. It may seem that we are back to 
the (historical) pre-inflationary times, where the form of the initial power spectrum of the 
Harrison-Zel'dovich type had to be postulated in order to reproduce the observations of the 
large-scale structures. This is somehow a very standard approach in physics in which one 
learns about the initial conditions of a system by observing its evolution 

Inflationary theories allowed to actually predict this spectrum, which makes them very 
predictive. We realize with this study that these inflationary predictions are very sensitive 
to the existence of an (classical) initial shear. To recover such a predictive power, we have 



to hope that a theory handling the dynamics of the universe on this scales jH, 0, 50 



or 



allowing to generate the shear [51[ will also provide a better understanding of the initial 
conditions. 



Indeed, one may wonder whether the arbitrary (pre- WKB era) conditions can be amplified 
and compete in amplitude with the ones of quantum origin seeded during the WKB regime. 
To estimate this, note that at early time (apart from the particular case a = '/r/2), 
behaves as —z"/zs and ul behaves as —z'l/zx- The expansion (]B2p then leads to the 
conclusion that, at early times. 
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whereas at late time, we deduce from our previous analysis that 



Thus, the solutions of Eq. fl3.23p has an oscillatory behaviour until the time when uj^ = Q 
and the arbitrary initial conditions are not amplified. For the tensor modes the situation is 
different because for some configurations of k, we can have uj\ ~ — at early time, and this 
leads to an exponential growth. However, since |u;a| "^C /rj with C = 0{1), this exponential 
growth is typically at most of order 



It follows that the arbitrary initial conditions either are not amplified or do actually grow, 
but in the latter case they are amplified by no more than a factor of order unity. We 
thus expect the transitional tachyonic behaviour and our ignorance of the initial state of 
the perturbations before the WKB regime not to significantly alter the validity of our 
prescription for the initial conditions of quantum origin set during the WKB regime. 

From a practical point of view, we can estimate the effect of the coupling functions K 
and D by changing their amplitude by hand. We have checked that they hardly affect 
the predictions that are presented in the following section. In particular, that teaches 
us that the main directional dependence of the power spectra is induced mainly by the 
directional dependence of the comoving wavenumbers (see Fig. [8]). We have also checked 
that when k/k^ei increases our prediction is similar to the one obtained in standard inflation. 

In conclusion, we can trust our prescription for modes larger than fcrcf and we have to 
assume (somehow as an observational input) that fcref < ^o- This sets a tuning on the 
primordial shear that we cannot explain with the model at hand. In this regime, we have 
checked that it is a robust prescription that is not affected by the unknown preexisting per- 
turbations that are completely arbitrary and that fixes the properties of the long wavelength 
modes. 

V. PRIMORDIAL SPECTRA: NUMERICAL EXAMPLES AND PREDICTIONS 
A. Definition of the spectra 

To set the initial conditions as previously detailed, we need to first solve numerically the 
system fl3.23H3.24l) . Because of the couplings, each of the three variables, v and /xa, will have 
components along the three independent stochastic directions, even if they were initially 
independent. For instance. 




^k(^) 



v^{\^, ?7)e„(k) + i;+(k, ?7)e+(k) + Vy^ (k, vi)e^ (k) 



and 



AtkA(?7) = /iAt,(k,?7)ej,(k) + /iAA(k, ?7)eA(k) + /iA(i_A)(k, ?7)ei_A(k) 
and we deduce from the properties of the random variables that 



Mr/K,(r/)) = (K(k,r^)|2 + K(k,r^)|2 + K(k,r/)|2) 5^^\\.-\^) . 



(5.1) 
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The power spectrum of the curvature perturbation at the end of inflation (once the shear 
has decayed away) is thus 



27rJ 



.P7^(k) = ^(|^;.(k,r/)|2 + |^;^.(k,r/)p 



(5.2) 



It can be checked, as expected, that for super-Hubble modes, 71 is conserved once the shear is 
neghgible. We thus perform our numerical integration in a time interval long enough so that 
the universe has been isotropized and all the observable modes have become super-Hubble. 
The power spectra of the gravity waves are defined analogously by 



^A(k) 



(5.3) 



At the beginning of the radiation era, the background spacetime can be described by a 
Friedmann-Lemaitre solution and the primordial anisotropy is now encoded on the statistical 
properties of the perturbations on large scales. It is thus convenient to decompose the power 
spectra on spherical harmonics according to 



and 



P,(k) = Mk) 



£=oo m=+l 

1=1 m=-e 



£=oo m=+i 



1 + E E rUk)YiM 



-1 m= 



(5.4) 



(5.5) 



The three functions fnik) and fxik) represent the power spectra averaged over the spatial 
directions, 

^-("^^ 47 • 



fx{k) 



The three series, rimik) and r^^(A;), characterise the deviation from statistical isotropy. 
Indeed, we expect the anisotropy to be negligible on small scales, that is 

rem{k) when k ^ k^ei , 

and that 

i^£m{k) — * when £ ^ 1 , 

the same being true to r^„^{k). Additionaly, because of the symmetries of the spectrum, the 
only non-vanishing coefficients are obtained for even i and even m. It can also be checked 
that these coefficients are real and that their values do not depend on the sign of m. Thus, 
we conclude that there are only 1 + i/2 independent real coefficients, 

e M i = 2i' ,m = 2m , m =0..J' . 

The gravity waves and curvature perturbation are also correlated so that 

{vk{r])fJ'Xk'{v)) = {\vvfJ'Xv\ + \vxfJ'Xx\ + \v{i-x)fJ'X{i-x)\) 5^^^(k- k') . (5.6) 
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B. Predictions 



To illustrate the signatures of a Bianchi / inflationary era, we consider a Bianchi spacetime 
with a = 7r/4. We fix the initial value Lfo/Mp = 3.25\/8tt ^ 16.3 for the inflaton field. 
This implies that the number of e-folds of the accelerating phase is — 67. We also 

set the reference wavenumber to k^^i — 157m ^ 157 10~^Mp. This corresponds to the 
definition (14.121) evaluated with the approximate cosmological constant solution fl2.33l) and 
(I2.38p . with the expression (12.571) for assessing r^:. 

We solve numerically the dynamics of the background and the evolution of the perturba- 
tions, as detailed in Appendix A, in order to compute the spectra defined in § IV A[ In order 
to understand their behaviour, we present: 

• Fig. [inj the evolution of the functions and fx{k) as a function of k for modes 
ranging from 3fcref to lOOfcref- This shows the evolution of the isotropic part, which 
dominates the small scales. 

We conclude that the curvature perturbation power spectrum has a spectral index 
ng — l — —0.032. For this model, 5 -C e and the modes depicted become super-Hubble 
when e is almost constant and e ~ 0.008 (see Fig. [6]). The expected spectral index 
in standard isotropic inflation is thus of order Ug — 1 = 25 — 4e ~ —0.032, in full 
agreement with our numerical computation. 

For tensor modes, we clearly see that P+ and Px differ on large scales and converge to 
the spectrum predicted by standard inflation on small scales. In particular, it can be 
checked that ut = — 2e ~ —0.016, in full agreement with our numerical computation. 

• Fig. [T71 a Mollweide projection of 'P7^(k) for different wavenumbers and for the scalar 
modes. This provides a visual intuition of the isotropization on small scales. 

• Fig. [T8j rgrnik) and r}^{k) as a function of k for the lowest multipoles (£ = 2). It 
quantifies the isotropisation on small scales, as was observed on the previous figure. 

Such predictions can be generated for any Bianchi / spacetime and up to an arbitrary 
multipole. In Appendix C, Ave provide another example, namely of the particular case 
a = tt/2 considered in Ref. |24l |. 

VI. CONCLUSION 

In this article, we have worked out the predictions of an anisotropic inflationary era for a 
generic Bianchi I spacetime. We have discussed the isotropization both at the background 
level and at the linear order in perturbation theory, both for scalar modes and gravity waves. 

Generically these spacetimes al way s enjoy a bouncing direction, apart from the particular 
case a = 71/2 considered in Ref. [2J]. (Note that the predictions for this case are in fact 
singular among the predictions since they do not converge uniformly when a — '/r/2). 

Since at early time, the modes are not in a WKB regime, we had to extend the standard 
procedure to fix the initial conditions (see § IIV Cp . We showed that the modes larger than 
/Cref always enter a WKB regime before they become super-Hubble during inflation. As we 
discussed, our procedure reproduces the standard one on small scales. 
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FIG. 16: Evolution of log [/ij (A:)] (left) and log[/A(A;)] for the two polarizations(right) as a function 
of log[A;//cref] for ipQ = 16.3Mp. The FL case is in dashed line. We also depict (in the inner right 
figure) the relative difference between the two polarisations, which shows that on small scales 
we recover that Px = -P+, as expected when isotropy is restored. On the upper line the initial 
conditions where fixed at tj(k) whereas in the bottom line they where fixed at r*. We see that 
for log {k/k^c{) > 0.8 the spectra are identical. For smaller wavenumbers, the WKB regime is too 
short to unambiguously fix the initial conditions. 



In the particular case where we tune the initial shear so that these initial conditions can 
be set unambiguously while still having an imprint of the CMB anisotropy, we presented 
the imprint of the primordial anisotropy in the power spectra of the gravity waves and 
curvature perturbation at the end of inflation. Two examples were studied but we can 
provide predictions for any Bianchi / universe. Note that these predictions were drawn 
by assuming that the slow-roll attractors were reached before the time the modes of 
observational relevance had exited the horizon. In the case of inflation with small number 
of e-folds before that time, it is not clear that this is a realistic hypothesis. If so, one would 
have to make the predictions in terms of trajectories, that is predictions that will depend 
on the initial conditions of the scalar field. This is not specific to Bianchi universes but to 



all models in which the inflationary period is short (40 



Concerning the initial conditions, two problems arise (see discussion in § lIVPp . First, 
there exists an early shear-dominated phase where the WKB approximation is violated. 
This forbids us to set the initial conditions as in a Friedmann-Lemaitre universe. As we 
showed, the sub-Hubble modes at the onset of the accelerating phase can be quantized 
because quantum fluctuations act at all times and can always source oscillatory solutions, 
when they exist. On the other hand, there always exist non-oscillating modes. These modes 
are expected not to be much smaller than the Hubble radius today (since there is no trivial 
imprint of the anisotropy on the CMB). This implies that there exists a cut-off scale above 
which we cannot predict the spectra from first principles, at least in the theoretical set-up 
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formalism has been designed when no interaction-free regime can be exhibited. 

In our analysis we have assumed the validity of general relativity up to the singularity. 
Indeed, we do not take this model for more than what it actually is, e.g. we cannot 
extrapolate it beyond the Planck or the string scale. There, more degrees of freedom 
have to be included and can change the dynamics. This could introduce an early chaotic 
phase jisl since Bianchi / models coupled to p-form fields have never ending oscillatory 
behaviour exhibited by generic string theory. Examples of an early dynamics have also 
been given in terms of Kalb-Ramond axion fields j49|, non-commutative geometry j5l| and 
recently loop quantum gravity [13] . Any of these developments may give a description of 
both the early phase and a procedure to fix the initial conditions. 



Coming back to inflation, our study demonstrates to which extent its predictions are 
sensitive to initial (classical) large scale anisotropies and that in the presence of a non- 
vanishing shear, it is impossible to define a Bunch-Davies vacuum in the standard way (see 



also Ref. [53| for a discussion of this issue in the standard picture and the arbitrariness on 



the choice of the initial state and its influence of the prediction of inflation). Indeed, if we 
tune the initial conditions such that the number of e-folds is large, none of the problems 
we address here will affect the observable modes, simply because only modes such that 
/c/fcref ^ 1 are observable and we have shown that in this limit we recover the standard 
inflationary predictions. If this is the case, one would have no observational imprint of the 
primordial shear. But our analysis and conclusions may be of some relevance for inflationary 
model building in the framework of string theory if the feeling that no large field model (and 
thus no large number of e-folds) can be constructed persists, an issue far beyond the scope 
of this article. 

Our analysis also shows the importance (and peculiarity) of the Friedmann-Lemaitre 
background in our theoretical predictions and on the quantization procedure, and it gives 
as well an explicit construction of the difficulties encountered when these symmetries do not 
exist. We showed that if the number of e-folds is large the inflationary predictions converge 
toward the isotropic predictions, hence demonstrating that they are robust in that regime. 
In the case of a small number of e-folds, they are very sensitive to the initial shear but also 
the theoretical construction is less under control. 

If the indication of the breakdown of statistical isotropy from the CMB were to be con- 
firmed and related to such an early anisotropic phase, then a new coincidence will appear 
in the cosmological models since one would need to understand why the characteristic scale 
is of the order of the Hubble radius today, k^cf ~ ko. 

From a more pragmatic attitude, the present work allows us to draw the CMB signatures 
of such an anisotropic early phase, for any Bianchi / universe. We stress that the general 
form of the initial power spectra are more general than those heuristically considered in 



the analysis performed in Refs. [35|, l5J] and go beyond the one derived for the singular 



case a = 7r/2 studied in Ref. [2J|. The existence of correlation between gravity waves and 
curvature perturbation, as well as the fact that the two gravity wave polarisations do not 
share the same spectrum, may also lead to specific signatures to be investigated. 
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APPENDIX A: INTEGRATING THE PERTURBATIONS 



We are interested in both the evolution of a cosmological mode and the decomposition 
of the shear according to this mode until it reaches the Friedmann stage. The covector of 
components fcj is constant by definition [see the discussion below Eq. fl3.ip ]. but in the Bianchi 
regime the vector /c* = Y'^kj is not constant since 'jij (and thus 7*-') is time dependent. This 
is why we used these covectors to label a mode. 

Once the background equations are solved for the functions Pi{t) and thus ai(t), jij is 
completely determined so that 

= E ^ (Al) 



and the associated normal vector is 



k- " ■ k- 



1. Evolution of the shear components 

The decomposition of the shear with respect to k also requires the construction of an 
orthonormal base |ei(t), e2(t), k(t)| that shall satisfy [2^ 

e\{t).e2{t) = ei(t).e'2(t) , (A3) 

where the scalar product is meant in terms of the spatial metric 7jj. 

In principle, it is only necessary to determine this basis at a given initial time, t\.a\t say, 
in order to extract the initial value of the shear components through Eqs. (13.8113. 91) 

5^(iimt) = (o'||,0'v1''^v25C^tX,0-t+)- (A4) 

Then, to determine the value of S at any time, we use that Eq. (12.211) implies (see Ref. [2^ 
for details) 

a;+27^a„=-2 5^a^„, (A5) 

a 

^va + 2^^va = 2^va^l| - E ^vb^xA-^ab , (A6) 

+ 2na^^ = 2 -^afe^va^vb > (A7) 
a,b 



where the matrix A^^^ is defined by 



which is manifestly symmetric in ah. It is explicitly given by 
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One interesting aspect of the system (lA5tlA7p is that it does not depend exphcitly on fcj. 
Once this system has been solved, we can deduce easily the functions (0;^,, cja, ^, ^) that 
enters the equations of evolution of the pertubative modes. 

In practice however, in the regime where a ^ 7Y, small numerical oscillations in the 
components of the shear are converted into huge numerical instabilities in the system (lA5t - 
IA7p . In order to avoid these numerical instabilities we turn to another method which we 
now describe. 



2. Systematic construction 



R {9/N) = 1 + J ■ fi— = 1 + Y^Jjfii — 



This method relies on the fact that the solutions of the system flA5HA7p are, in a general 
sense, given by 

cr|| = aijky , a^a = CTijk'el , a^x = o^ije'^ (AlO) 

as a function of time. 

Instead of solving flA5tiA7p after having extracted the components of the shear at a fixed 
initial time, we can determine the dynamics of the shear components through the dynamics 
of the basis vectors, the polarisation tensor and the shear. Once the are determined 
numerically, the shear is known. 

In order to determine the base vectors {ei(t), e2(t)}, we first start from a triad 
{ex(t),ey(t),ez(t)} aligned with the xyz-piopei axis. Then we introduce three Euler an- 
gles to rotate this triad to any given direction. 

Let us recall how to deal with rotations in spaces endowed with a general (non orthonor- 
mal) metric. Any rotation matrix about a unit vector n can be constructed by means of 
infinitesimal rotations of the form 

e 

'N ' ' ' "'"'iV 

where 6/N is some infinitesimal angle. The Jj are the generators of the group of rotations. 
In particular, if we want to specify a rotation about the 2;-axis, then the unit vector aligned 
with this axis has components 

n,= (0,0, e'^^) , 

and we have 

/ n\N 

(6) = Jhn M + e-^'J^-j = exp [Oe-^U,] . 

It can be shown that in order to conserve the scalar product, the generators must satisfy 
{Jk)ij = -^kij, where 

^ijk is totally antisymmetric normalised such that €123 — v/det(7i,) = 1. 
A rotation matrix around the 2;-axis then requires to use {JzYj = {lY^{Jz)kji which is 
explicitly given by 

/ -e-2ft 0\ 
{J.)\ = . (All) 

\ / 

Consequently, taking into account the constraint ( 12. 5p . i.e. Ylil^i — O5 rotation about 
the 2;-axis can be written as 

cos(^) -e(^2-/3i) sin(^) g 

[Rz{0)][,= \ e^f^^-f^^^ sm{9) cos(^) | . (A12) 
1 
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Similarly, for a rotation about the y-axis, we use the generator 

/ e-2ft \ 

{Jy)\ =\ (A13) 

\-e-2/33 / 

to get 

/ cos(^) e(^3-^i) sin(^) 

[Rym% = 1 I . (A14) 

y _e(/3i-/33)sin(^) cos(^) 

Explicitly, the components of the triad from which we start are 

(e.r=|^'oj, (e,r ^ 1^ e ^/^^ j , (e,)^ ^ | |. (A15) 

The three Euler angles 7) are then used to rotate this triad according to, respec- 

tively, the X, y and 2;-axis. Explicitly, the triad after rotation is given by 

{e,r ^ RM',Ry{f3)\R.{ay^{e,r, 

{e,r ^ RMyRy{P)\RM,{eyY. 

{ukY ^ RM'^RyW)\RMpi^zY. (A16) 
which determines the orthonormal basis {ei,e2,n/c}. This prescription is still incomplete 
since our problem also requires the covector Uj. to be equal to the covector k during the 
whole evolution of the system, i.e. 

{uk)i{t) = kit) , Vt . 

Since ki should not depend on time, we need to determine the Euler angles as a function of 
time such that, for any two times t and t', we have [uk{t)]i = f(t,t')[uk(t')]i. This condition 
is satisfied provided 

tan(7) = tan(7^) exp [(A - /^s)] (A17) 



tan(/3) = exp 



cos(7) 



(A18) 



where 7/ and (3f are the angles which give the final direction that we consider when the 
shear has vanished. 

Additionally, our problem also requires the condition 

[{e2y]' [{e,U = [{e,y]' [{e^U , 

as a prescription for the choice polarisation basis to be continuous. It is fulfilled if 

a' = - cos{p)-f' . (A19) 

This equation can be integrated with the help of Eq. (1A17|) and Eq. (lAlSp . Consequently, 
provided we know (/?/,7/) describing the final orientation of a mode, we can determine 
{a(t),^{t), /3{t)) which will select the triad adapted to the mode at any time. 

In conclusion, the set of equations (lA16tlA19l) gives a complete description of the time- 
evolving triad necessary to extract the components of the shear according to Eq. flAlOP at 
any time. This is the procedure we have implemented in our numerical calculations, and is 
robust to numerical errors. 
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APPENDIX B: SLOW-ROLL EXPRESSIONS 



We give here, for the sake of completeness, the exphcit expressions for ^sV-^s, ^a/^a, 
and n in terms of the slow-roll parameters. Setting 



_ ^11 _ ^tA 



we have 



^A 



W{2 + e 



5 - 45 + 2e 



(3 - 5f 
3-e 



[-3 - 65 + e (7 + 45 - 4e)] + 2x^ (3 - ef 
+ (y^) [6e + - 4e5 + {-Ae^ + 4e5) + 2x^ (e^ - 6e)] | 



<^ 2 - e + 6w; - 2e 



w + 'i 



yl 



1 — w 



(e - 3) x^ 



1 + 2u; + 6 



yl 



1 — w 

2/A 



+ 12Z/(l-A) + 1 

w' 6 



1 — U7 



yl 



n\i-w 



en' 



n\i-w 



6(3-e + (3-e)x2) f | 



y+y> 



w 



(Bl) 



(B2) 



(B3) 
(B4) 
(B5) 



APPENDIX C: PARTICULAR CASE a = Tr/2 

As we stressed, the case a = 7r/2 is particular in the sense that this is the only Bianchi / 
models which does not have any bouncing direction. 

Since this model enjoys a planar symmetry (conveniently chosen here to be in the xy- 
plane), the component cr^x of the shear vanishes for any choice of the Euler angles provided 
that a = 0, which corresponds to the symmetry around the 2;-axis. This implies that 
oc (Tx also vanishes at all times. The function H behaves as cr+cix and for the same 
reason also vanishes. In conclusion, when (and only when) a = -7r/2, we have 







n = o . 



This is the reason why it was possible to construct a "simplified" perturbation theory in 
Ref. 2J] in splitting in 2 + 1 degrees of freedom from the start of the analysis. As we saw, 
this is not generic. 

Figure [12] now shows that 

I- 
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at early time which imphes that the modes were in fact decoupled. Figs. [20] and [21] consider 
the validity of the WKB condition. In particular, one can show that on sub-Hubble scales 

7" h"^ 7" 

fs ^ ^ _£ 

Zs aj Zx aj 

We have checked that there exists only one time redefinition which leads to adiabatic canon- 
ical variables and that agree with all the conclusions in Ref. j23 |. 

Following our quantization procedure, we compute the predictions for the power spectra 
(Figs. [22] to [24|) . Note that the planar symmetry is obvious on the Mollweide projection, 
that shall be compared to the generic case presented on Fig. [TT] This illustrates how a 
detection of anisotropy may permit us to reconstruct the particular Bianchi / structure (if 
any !) during inflation. 

In this particular case, there always exist a WKB regime We can thus define our initial 
conditions in the standard way. Following the procedure described in this article, we can 
compute the power spectra of this solution. We present here the same figures as in § IV B[ 
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FIG. 19: Evolution of H^/W^ for the three orthogonal modes aligned with the and z-axis 
(represented by three different colors). We have considered here a generic Bianchi spacetime with 
a = 7r/2. 
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FIG. 20: Evolution of IQ^^^/cj^l for three different modes, each of them ahgned with one of the 
three ortoghonal directions and with the same modulus lO/crcf at the end of inflation. We have 
considered a generic Bianchi spacetime with a = Tr/2. 
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FIG. 21: Evolution of IQjf^^ /^xl for A = + (left) and A = x (right) for various modes with the 
same modulus 10A:rcf at the end of inflation. We have considered a generic Bianchi spacetime with 
a = 7r/2. 
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FIG. 22: Evolution of fR{k) (left) and fx{k) (right) as a function of log[A;/A;i.cf]. The FL case is in 
dashed line. 
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FIG. 23: Mollweide projection of the ration between V-jiili) and its value in the FL case expressed 
in percentage for log[k/k^cf] = 1/2, 3/4, 1, 3/2 from left to right, top to bottom. The x-axis as 
been choosen as the vertical axis 
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FIG. 24: Evolution of rim{k) (left) and r^mi^) ('^ight) as a function of log[A;/fcref] for the lowest 
multipoles (1 = 2). The x-axis has been used to perform the spherical harmonics decomposition. 
As a result of symmetries, only the coefficients with m = do not vanish. On the right the red is 
the + polarisation and the blue is the x polarization. 
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